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Microlensing as a probe of the Galactic structure 
20 years of microlensing optical depth studies 
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Abstract Microlensing is now a very popular observational astronomical tech- 
nique. The investigations accessible through this effect range from the dark 
matter problem to the search for extra-solar planets. In this review, the tech- 
niques to search for microlensing effects and to determine optical depths 
through the monitoring of large samples of stars will be described. The conse- 
quences of the published results on the knowledge of the Milky- Way structure 
and its dark matter component will be discussed. The difficulties and limita- 
tions of the ongoing programs and the perspectives of the microlensing optical 
depth technique as a probe of the Galaxy structure will also be detailed. 
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1 Introduction 

In 1936 A. Einstein published a half-page note in Science entitled "lens-like 
action of a star by the deviation of light in the gravitational field" |48j . with 
this hnal comment: "There is not great chance of observing this phenomenon. " 
According to J. Renn [93] the effect was already predicted by A. Einstein as 
early as 1912, before completion of the general theory of relativity. S. Liebes 
[71] considered again the properties of the gravitational lenses in 1964, and 
studied the possibility to detect high magnification events. Furthermore, this 
author already mentioned the possibility to detect invisible compact objects 
through microlensing. In 1986, B. Paczyhski [5T] pointed out the possibility of 
using the gravitational microlensing effect to detect massive compact objects 
of the Galactic halo in the direction of the Magellanic Clouds. 
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In 1989, the raising question amongst astrophysicists was the nature of 
dark matter within the Milky- Way halo. WIMPSLj search experiments were 
already looking for exotic dark matter, but it was realized at this epoch that 
the question of today's nature of the majority of the baryons (at z = 0) — still 
pending — could be addressed through the microlensing technique. From 1989, 
several groups started survey programs to search for compact halo objects 
within the Galactic halo. The challenge for the EROS and MACHO teams was 
to clarify the status of the missing hadrons in our own Galaxy. In September 
1993, three teams, EROS [SJ, MACHO and OGLE QTT] discovered the first 
microlensing events in the directions of the Large Magellanic Cloud and the 
Galactic Center. Since these first discoveries, thousands of microlensing effects 
have been detected in the direction of the Galactic Center (GC) together with a 
handful of events towards the Galactic Spiral Arms (GSA) and the Magellanic 
Clouds. 

Probing Galactic structure with microlensing: Microlensing has proven 
to be a powerful probe of the Milky- Way structure, and not only for the hid- 
den baryonic matter. Searches for microlensing towards the Magellanic Clouds 
(LMC, SMC) and M31 provide optical depths through the Galactic halo, al- 
lowing one to study dark matter under the form of massive compact objects. 
Searches towards the Galactic plane (Galactic Center, GC and Galactic Spiral 
Arms, GSA) allow one to measure the optical depth due to ordinary stars in 
the Galactic disk and bulge. Kinematical models as well as mass functions can 
also be tested through the event duration distributions. In a few special cases 
called "exotic" events, more specific information on the lensing system has 
been obtained; unfortunately they don't yet provide sufficiently large samples 
to extract reliable information on the Galactic structures. 

This article reviews the optical depth measurements, which are the instan- 
taneous probabilities for a point-source of a target to be magnified by a factor 
larger than 1.34. The description of the observational devices and of the analy- 
sis techniques will be given, followed by a discussion about the main difficulties 
and limitations which can be overcome in future surveys. The consequences 
of the optical depth measurements on the Galactic structure knowledge will 
also be reviewed, with the first connections that can be established between 
the kinematics of the Galactic structures and the duration distributions. We 
will conclude with a prospective on the future surveys and a description of the 
potential of microlensing searches in infrared and of "fine" searches, sensitive 
to exotic events, providing valuable extra information. 

2 Microlensing basics 

Gravitational microlensing effect occurs when a massive compact object passes 
close to the line of sight of a star and produces gravitational images that 
are not intercepted by the massive object (no eclipse), but that cannot be 
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Fig. 1 Principle of the microlensing effect: As the deflector D of mass M moves with a 
transverse relative speed Vt, the impact parameter u(t) changes with time, and so does the 
magnification of the source. 



separated in telescopes (no multiple images) (see Fig. [IJ . The detection of such 
a coincidence is then possible only if a measurable magnification variation 
occurs during the observation time. The details on the formalism associated 



with this phenomenon can be found in [100 and we will only provide here a 
few basic tools. In the approximation of a single point-like object acting as 
a deflector on a single point-like source, the total magnification of the source 
luminosity at a given time t is the sum of the contributions of two images, 
given by 

Ait) = < t? + 2 , (1) 

u{t)^u{ty + A 

where u(t) is the distance of the deflecting object to the undcflcctcd line of 
sight, expressed in units of the "Einstein Radius" Re- 
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Here G is the Newtonian gravitational constant, D$ is the distance of the 
observer to the source and xD$ = Dl is its distance to the deflector of mass 
M. The motion of the deflector relative to the line of sight makes the mag- 
nification vary with time. Assuming a deflector moving at a constant relative 
transverse speed vt, reaching its minimum distance uq (impact parameter) to 
the undeflected line of sight at time to, u(t) is given by 

u{t) = ^u 2 + {(t-h)/t E ) 2 , (3) 

where t E = R e /vt, the "lensing time scale", is the only measurable parameter 
bringing useful information on the lens configuration in the approximation of 
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simple microlensing: 
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2.1 Microlensing event characteristics 



The so-called simple microlensing effect (point-like source and point-like lens 
with uniform relative motion with respect to the line of sight) has some char- 
acteristic features which allow one to discriminate it from any known intrinsic 
stellar variability : 

— The event is singular in the history of the source (as well as of the deflector). 

— The gravitational origin of the effect implies that the magnification is in- 
dependent of the color. 

— The magnification is a known function of time, depending on only 3 pa- 
rameters (u m j n , to, £e), with a symmetrical shape. 

— As the geometric configuration of the source-deflector system is random, 
the impact parameters of the events must be uniformly distributed. 

— The passive role of the lensed stars implies that their population should 
be representative of the monitored sample, particularly with respect to the 
observed color and magnitude distributions. 

This simple microlensing description can be broken in many different ways: 
the lens may be a double system [75] . or the source may be an extended object 
[120j , or the relative motion with respect to the line of sight may deviate from 
a uniform motion due either to the rotation of the Earth around the Sun 
(parallax effect) [55] [55], or to an orbital motion of the source or of the lens 
around the center-of-mass of a multiple system [75] . All these deviations have 
been observed and will be discussed further in this issue (see the contribution 
of M. Dominik [IS]). The interest of these "exotics" comes from the extra- 
information on the lensing configuration that can be obtained in such specific 
cases. 



2.2 The observables: optical depth, event rate and ts, distribution 

The optical depth up to a given source distance Ds is defined as the instanta- 
neous probability for the line of sight of a target source to intercept a deflec- 
tor's Einstein disk, that corresponds to a magnification A > 1.34. Assuming 
that the distribution of the deflector masses is described by a density function 
p(Dl) and a normalized mass function dni(Z?L, M)/dM, this probability is: 

where Be = Rr/Dl is the angular Einstein radius of a lens of mass M located 
at Dl. The second term of the integral is the differential number of these 
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lenses per mass unit and per solid angle. As the solid angle of the Einstein 
disk is proportional to the deflectors' mass M, this probability is found to be 
independent of the deflectors' mass function: 

t{D s ) = / x(l - x)p(x)dx , (6) 



where p{x) is the mass density of deflectors located at a distance xD$- The 
mean optical depth towards a given population defined by a distance distri- 
bution dns(Ds)/dDs of target stars is defined as 

The measured optical depth associated to microlensing events observed in 
a population of N b s stars monitored for a duration AT b s is given by 

T N obs AT obs 2 ^ e(t E ) ' 1 ' 

events y ' 

where e(t^) is the average detection efficiency of microlensing events with a 
time scale (see the precise definition in Sect. 14. 4| ). 

The calculated optical depth does not depend on the deflectors' mass func- 
tion, but the measured optical depth, that takes into account the mean de- 
tection efficiency e^e), can be biased by this mass function, in particular 
because the detection function vanishes for very short or very long duration 
events. This fact makes impossible a perfect compensation of the inefficiencies. 
If many deflectors are light (resp. heavy) enough to produce extremely short 
(resp. long) duration events that cannot be detected, the measured optical 
depth will clearly be underestimated. This is why collaborations often indi- 
cate that their results are valid within a given duration domain. The control 
of the detection efficiency — that is needed to get reliable optical depths and 
rate estimates — is one of the most delicate aspects of the survey searches. 

Event rate: Contrary to the optical depth, the microlensing event rate de- 
pends on the deflectors' mass distribution as well as on the velocity and spatial 
distributions. The global measured event rate, corrected for the detection effi- 
ciency, is 

r = ^ZAfZ x JC K^y (9) 

This event rate and the duration distributions are used to constrain the mass 
function and the kinematics of the lensing structures [53] [121] . Such stud- 
ies will only be briefly mentioned in this review, as the available efficiency 
controlled statistics are still very limited for quantitative comparisons. 



2 The product N ), s X AT oba is called the exposure. 
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3 Observations 

Three types of programs can be distinguished. The two first ones, that can be 
called "survey programs" , monitor large samples of stars that are either iden- 
tified in the images (catalog survey programs) or that are not identified (pixel 
survey programs). The third type can be considered as "follow-up programs" 
that finely monitor individual events detected by the early warning systems 
of the survey programs. In this review, we will only discuss the catalog survey 
programs with brief mentions of the other types that are discussed in details 
in the other contributions of this volume; in particular, results from the M31 
pixel surveys are discussed in the Calchi Novati's review [35], and the follow-up 
programs are reviewed by M. Dominik |45) . 

Table Q] lists the characteristics of the devices used in the catalog survey 
programs. The common requirements of these surveys are the monitoring of 
a very large number of stars, with the best time-sampling and the best pho- 
tometric resolution that can be achieved. These requirements drive the main 
specifications of the instruments: they all monitor wide crowded fields with 
broad passband filters -to get a large photon flux. The achromaticity of the 
microlensing events further supports this use of wide passband filters; since 
the ratio of luminosities does not depend on the wavelength, the signal sig- 
nificance increases with the photon flux. These surveys use in general at least 
two different passbands to make easier the identification of achromatic events. 
Figures [3] and [5] show the fields monitored by the different teams towards the 
Galactic plane and the Large Magellanic Cloud. 

EROS: Experience de Recherche d'Objets Sombres. After a first phase 
started in 1990, using ESO-Schmidt (lm) telescope photographic plates dig- 
itized with the MAMA machine and a 16 CCDs camera [25] [23], the EROS 
team used the lm (F/5) MARLY telescope installed at La Silla observatory 
(Chile) from July 1996 to February 2003, with a dichroic beam-splitter allowing 
the simultaneous imaging in two wide pass-bands (EROS-bluc and EROS-red) 
of a 0.7° (a) x 1.4° (5) field [25]. Photons were collected by two cameras equipped 
with eight 2K x 2K LORAL CCDs each. Exposure times ranged typically from 
2 to 5 minutes and were optimized to maximize the global significance of the 
photometric measurements taken during microlensing magnifications. EROS 
monitored LMC and SMC, the Galactic Center (GC), and four directions in 
the Galactic Spiral Arms (GSA). 

MACHO: MAssive Compact Halo Objects. From June 1992 to January 
2000 the MACHO team used a 1.27 m (F/3.9) dedicated telescope [67], also 
equipped with a dichroic beam-splitter and two cameras of four 2K x 2K LO- 
RAL CCDs each [iTj2] [75]. This setup, installed at Mount Stromlo (Australia), 
provided simultaneous blue and red band photometry and covered a field of 
0.7° (a) x 0.7° (S) with a pixel size of 0.62 arcsec. The LMC, SMC as well as the 
Galactic Center have been monitored for microlensing searches. The observing 
strategy changed during the project and the sampling of the Galactic fields 
varied from a few observations per year to daily observations. 
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Table 1 Experimental setups devoted to microlensing surveys, with the telescope diameter, 
the list of the filters, the total number of pixels, the field covered by the detector, the list 
of surveyed targets with the corresponding monitored fields, numbers of monitored stars, 
approximative average sampling and observation epochs. The indicative sampling values are 
averaged over the observation seasons or during the observation nights. 

Note: the MACHO SMC data were lost during the Mount Stromlo fire of 18 January 2003. 



Set-up 


telesc. 
filters 


pixels 
xlO 6 


detector 
field (°) 


target 


field 
deg 2 


stars 
XlO 6 


approx. 
sampling 


obs. 
dates 


EROS1 
plates 


1 m 
R/B 


1000 


5.2x5.2 


LMC 


27. 


4.2 


3 days 


90-93 


EROS1 
CCD 


0.4 m 

R/B 


4 


0.4x1.1 


LMC 

SMC 


0.44 
0.44 


0.1 
0.1 


20 min. 
20 min. 


91-94 
93-95 


EROS 2 


1 m 
I/V 


2x32 


0.7x1.4 


LMC 
SMC 
GC 
GSA 


84. 
9. 
66. 
20.1 


29.2 
4.2 
60. 
12.9 


3 days 
3 days 
3 days 
3 days 


96-02 
96-02 
96-02 
96-02 


MACHO 


1.27 m 
R/B 


2x16 


0.7x0.7 


LMC 
SMC 
GC 


14.7 
3. 
48. 


11.9 
2.2 
50.2 


4 days 
3 days 
3 days 


92- 99 

93- 96 
92-99 


OGLE I 


1 m 
I 


4 


0.25x0.25 


GC 


1.25 


1.5 


3 days 


92-95 


OGLE II 


1.3 m 
I/B/V 


4 


0.24x0.95 
drift-scan 


GC 
LMC 
SMC 


11. 
4.7 
2.4 


20.5 
5.5 
2.2 


3 days 
3 days 


97-00 
96-00 


OGLE III 


1.3 m 

I/V 


64 


0.58x0.58 


LMC 
SMC 
GC 


38. 
13. 
35. 


total 
120. 


2 days 


01-09 


DUO 
plates 


1 m 

R/B 


1000 


5.2x5.2 


GC 


27. 


13. 


1 day 


94 


MOA1- 
caml 


0.61 m 
B/R 


9 


(l.Xl.)/4 


LMC 
SMC 


3.5 
2.8 


1. 
0.4 


3 hours 
3 hours 


95-98 
95-98 


MOA1- 
cam2 


0.61 m 
B/R 


24 


0.92x1.38 


GC 
LMC 
SMC 


20. 
20. 
10. 


5. 
3. 
1. 


2 hrs 
1 day 
1 day 


98-05 
98-05 
98-05 


MOA2- 
cam3 


1.8 m 
R/V/I 


80 


1.32x1.65 


GC 
LMC 
SMC 


48. 
31. 
4.4 


50. 
50. 
5. 


1 hr 
1 hr 
1 day 


05- 
05- 
05- 


Super- 
MACHO 


4 m 
VR 


64 


0.57x0.57 


LMC 


22. 


100. 


2 days 


01-05 



OGLE: Optical Gravitational Lensing Experiment. After a first phase 
started in f992 |110j . making use of a lm (F/7) telescope in Las Campanas 
(Chile) to monitor the Galactic Center, the OGLE team has built its own 
dedicated instrument, using a 1.3m (F/9.2) telescope, equipped with a thinned 
CCD camera of 2K x 2K pixels and a filter wheel with standard UBVRI filters 
|108j . The system operated in drift-scan mode with the telescope drifting in 
declination at a rate of a few arc-seconds per time second. In these conditions 
a single image covered a field of 0.25°(a) x 0.9° (8) with a pixel of 0.4 arcsec. 
LMC, SMC and the Galactic bulge have been monitored with this new setup 
from 1997 to 2000. Each field was observed every third night in I-band by 
OGLE-II, and every 11-th night in the V-band. OGLE-III had a 8K x 8K 
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pixel camera, made of 8 CCDs. Observations were performed towards SMC, 
LMC and the Galactic bulge Q2Z] from June 2001 to May 2009. OGLE-IV is 
made of 32 E2V CCDs with 2K x AK pixels each and received its first light 
in September 2009. 

DUO: Dark Unseen Objects. This search for microlensing towards the 
Galactic Center, also using the ESO-Schmidt telescope photographic plates 
digitized with the MAMA, was performed in 1994 [BJ. 

MOA: Microlensing Observations in Astronomy. This team started mi- 
crolensing searches with a 0.6m telescope located at mount John Observatory 
(New-Zealand), equipped with an array of 3 x 3 non-buttable IK x IK CCDs 
[50] . Two different focal lengths were used during this first stage. The second 
camera had three 2K x AK pixels thinned SITc CCDs |119j . The team is now 
using a new 1.8m telescope, with ten 2K x AK pixels E2V CCDs [§5] and it 
performs a real-time difference imaging analysis as described in 33 . 

SuperMACHO. From October 2001 to 2005, this collaboration used the 
4m CTIO Blanco telescope during 150 half-nights, equipped with eight SITe 
2K x AK CCDs and a 250nm passband VR filter, to monitor all star types 
in the LMC |93) . Considering the serious supernovae contamination, comple- 
mentary observations triggered through an alert system were used to identify 
microlensing events. 

4 Analysis techniques 

The data analysis of the surveys are all based on the same principles : Firstly, 
light curves (i.e. flux versus time curves) for cataloged stars or for emerging 
pixel clusters are extracted from the raw data, using automatic astrometric 
and photometric alignment procedures; secondly, each light curve is subjected 
to a microlensing search algorithm based on the expected characteristics of 
the events. 

4.1 Photometry 

Simple aperture photometry is not adapted for crowded fields. The microlens- 
ing pioneers have used or developed Point Spread Function (PSF) photometry 
techniques to extract the luminosity of each star in crowded fields ([50] and 
references therein for the photographic plates, [21] for the CCD). A simul- 
taneous fit of the PSF is performed on stars that are not completely separated. 
The following generation of microlensing teams have developed Difference Im- 
age Analysis techniques (DIA) [7] |115j [70] [33] that are very demanding in 
computing ressources, and are now usable thanks to the advances in com- 
puters. The base of this technique is the construction of the image difference 
between a current and a reference images. A convolution kernel is calculated 
and used to convolve the reference image to match the seeing of the current 
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image. The two images are then subtracted and the variability are searched 
from the study of the residuals. This technique has the advantage to make de- 
tectable the magnification of objects that are not detected when unmagnified. 

4.2 Production of light-curves 

Catalogs of objects are extracted from reference images, produced by the co- 
addition of series of the best quality images. For PSF photometry, the light- 
curves are produced after geometric and photometric alignment of current 
images on reference images. For DIA technique, a list of variable objects is 
extracted from the subtracted images and subsequently monitored in the other 
images. 

4.3 Selection of microlensing candidates 

The general philosophy for this discriminant analysis is as follows: if the start- 
ing point is a series of light-curves, the first stage consists in selecting the 
less stable objects, then to search for microlensing candidates through the 
topological characteristics of the expected light-curve. All analysis include a 
cleaning procedure, where bad quality images (taken under poor atmospheric 
conditions, with guiding problems...) are removed, sources with problematic 
environment (such as dead pixels or bright neighboring stars) are discarded 
from the reference catalog, and the identified aberrant measurements are sup- 
pressed. Various algorithms of discriminant analysis have been developed to 
pre-select light-curves with single positive fluctuations, visible in each color (in 
the case of multi- wavelength observations). Microlensing fits are performed on 
the pre-selected light-curves and the final selection is usually based on the 
fitted parameters. 

The last selection stage is the background rejection and the residual back- 
ground estimate. This stage has long suffered from being somewhat underes- 
timated. Apart from the known experimental artifacts such as contamination 
by egrets or those induced by the SN1987A echoes |118j . a few categories of 
objects that could mimic a microlensing event have been identified by the 
search teams over the years; the microlensing science would certainly benefit 
from a better understanding of these fakes: 

- The so-called "blue bumpers" [T5], bright blue variable stars that show 
recurrent bumps. Amongst them, Be type stars that are known to be unstable 
(case of the EROS1-1 candidate [27] ) 

- Supernovae that explode in galaxies behind the monitored target [TT] . 

4.4 The efficiency control 

The optical depth calculations require the knowledge of the detection effi- 
ciency as a function of £e 5 averaged over the impact parameter space, over 
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the observation period and over the source apparent magnitudes and colors. 
This detection efficiency is defined as the ratio of the microlensing events that 
satisfy the selection requirements to the number of events with uq < 1 whose 
magnification reaches its maximum during the observing period. Its calcula- 
tion needs the simulation of events in a domain of the parameter space that 
exceeds by a large amount the expected domain of sensitivity (in EROS, uq 
was generated up to 2, 1 day < tg < 900 days, t was generated from 150 
days before the first observation to 150 days after the last one). The simplest 
technique to simulate microlensing-like light-curves consists in superimposing 
simulated events on measured light curves from an unbiased sub-sample of the 
monitored catalog. Events are simulated as point-source, point-lens, constant 
velocity microlensing events. This technique automatically takes into account 
the real time sequence, the source variabilities and the photometric systematic 
or accidental distortions that are preserved when combining with a simulated 
microlensing through rescaling procedures. The simulation of the microlensing 
in series of synthetic images (instead of light-curves) have also been made by 
[T5] and |105j , followed by the full photometric reduction and light-curve pro- 
duction. MOA used an intermediate technique of superimposing a microlensing 
effect on real subtracted images [104] . 

When the monitored population has a large distance dispersion, then the 
microlensing detection efficiency for a given stellar type depends on the dis- 
tance, since it depends on the apparent magnitude. This makes the efficiency 
to be ultimately correlated with the optical depth itself. As the microlensing 
detection efficiency essentially depends on the apparent magnitude and possi- 
bly on the color, the mean optical depth up to a given magnitude (or for a 
given catalog) r cat is the concept that has to be used in the case of microlens- 
ing of populations with a wide distance distribution such as those belonging 
to the Galactic Spiral Arms [57]- Such a measured optical depth r cat can be 
compared with the predictions derived from a lens and source distribution 
model using expression (0, where the observed source distance distribution 
dns(Ds)/dDs corresponds to the cataloged sources (taking into account the 
star detection efficiency). 

4-4-1 Blending 

Due to source confusion, single apparent cataloged objects can be made of sev- 
eral stars. This fact, called blending, is one of the major sources of uncertainty 
for the optical depth determination in crowded fields, as it affects not only the 
microlensing detection efficiency but also the effective number of monitored 
sources N Q b s that enters in expression © [H] [5] [116] [55] . When some blend- 
ing is taking place, the reconstructed stellar flux contains contributions from 
more than one star. For this reason, the effective number of monitored sources 
N b s that should enter in expression ([5]) is larger than the number of cataloged 
stars. Also, as only one of the components of the blend is magnified, this sit- 
uation makes the reconstructed apparent magnification of the blend smaller 
than the real magnification of the lensed component. As a consequence, the 
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detection efficiency appears smaller than for the unblended case. The impact of 
blending on the optical depth estimates has been studied by the observational 
teams using different approaches; as an example [3] [55] did a study based on 
synthetic images and on real images enriched with randomly located synthetic 
stars. In this type of simulation, it is usually assumed that there is no spa- 
tial correlation between the blend components and that the blend results only 
from a random coincidence. As the existence of correlation should increase 
the blending, this should be systematically checked and taken into account 
if necessary for the simulations. Other studies [T7] |117j [57] compare HST 
images with the corresponding monitored sub-fields. For all these studies, the 
critical ingredient is the underlying luminosity function. Fig. [5] shows better 
than in words what is the blending in a Galactic plane field, and illustrates the 
relation between the (HST) luminosity function and the content of an EROS 
object. The HST images that are deeper by at least 2-3 magnitudes than the 
ordinary survey catalogs and have a PSF of ~ 0.1" — 10 times smaller than 
the terrestrial images — allow one to extract relevant statistical data in the 
underlying ingredients of the objects cataloged by the surveys. 

A correction to the efficiency estimates is derived from these studies. There 
is some compensation between the loss of event- by- event efficiency and the gain 
in the effective number of monitored objects, that limits the impact on the op- 
tical depth estimates. Finally, the systematic uncertainty on the optical depth 
due to blending is generally estimated to be smaller than 10%. Decreasing this 
uncertainty needs a careful estimate of the blending corrections. The most 
simple way to minimize the impact of blending is to keep only microlensing 
of red giant (bright) stars. This option has been used by all the search teams 
towards the Galactic Center, and by EROS and OGLE towards the LMC. 

4-4-2 The efficiency uncertainties due to complex events 

The efficiency estimates are based on the simulation of simple microlensing 
events. It has been estimated that about 10% of the real microlensing events 
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have a more complex evolution such as a double lens configuration or a paral- 
lax effect [S3] . In the later case, the correction is usually small, but in the first 
(rare) case, the mean detection efficiency may be seriously affected by the dis- 
torted shape of the light-curve [53] . The impact of this detection modification 
is sometimes accounted for as an extra-uncertainty of ~ 10% on the detection 
efficiency. 



4.5 Complementary observations and statistical tests 

These a posteriori tests are crucial not only to check and improve the back- 
ground rejection, but also to gain confidence in the statistical distributions. 
The first microlensing candidates have been investigated in details through 
complementary photometric and spectroscopic measurements. A non negligi- 
ble fraction of the early candidates (with low signal to noise) has been rejected 
after publication because of the late discovery of anomalies such as repetition 
of bumps {e.g. the so-called "blue-bumpers" 12 ) or baseline variability |22) . 
In connection with optical depth studies, systematic statistical tests are per- 
formed on the populations of detected microlensing candidates to check that 
the impact parameter uo and the maximum magnification time to distributions 
— with uniform prior — are only biased by the detection efficiencies, and that 
this population is representative of the monitored population after correcting 
for the detection efficiencies. 



5 Some limitations and difficulties in the optical depth 
determination and future improvements 

5.1 Statistical limitations 

The impact of poor statistics is worsen by the dispersion of the microlens- 
ing event durations; since the variability of the durations induces a strong 
variability of the individual contributions to the optical depth, the statistical 
uncertainties on this parameter are indeed larger than expected from Pois- 
son statistics. This is taken into account by the search teams by using the 
procedure described in |64) . 

One of the statistical limitations in the microlensing surveys is the available 
number of sources that can be monitored. As an example, optical depth mea- 
surements towards obscured regions of the Galactic plane are very strongly 
limited by the visibility of the sources. The only way to improve the situa- 
tion is to make infrared surveys [SB]. The VVV (VISTA Variables in the Via 
Lactea) program [114] , one of the vistaI [113] large survey programs, aims at 
searching for variabilities towards the Galactic plane and should have a major 
impact for the microlensing searches towards the dusty regions. 



3 Visible and Infrared Survey Telescope for Astronomy 
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Another statistical limitation — more difficult to overcome — comes from 
the weakness of the signal itself; more exposure (N b s x AT Q b s ) is needed to 
compensate for the low microlensing rate towards the Magellanic Clouds. Su- 
perMACHO [93] which monitored the largest possible number of stars within 
the LMC was specifically devoted to this search. 



5.2 Information from simple events 

The ordinary microlensing events provide little physical information on the 
lens configuration. Any technique improving the photometric accuracy through 
the use of larger telescopes or by achieving high sampling rates increases the 
chance to get some extra-information by the detection of small deviations 
with respect to the simple microlensing light-curve. International networks of 
relatively small telescopes linked to the alert systems of surveys like PLANET0 
[85] MicroFUPN0 [77] or MINDSTErJl [78] have developed the second approach. 

The LSST0 project [72] will certainly dramatically change the landscape for 
the photometric accuracy, as the telescope has an equivalent diameter of 6.5m. 
Tens of thousands of events should be detected each year with a sampling rate 
of once every few nights. Considering the photometric accuracy, parallax and 
extended source effects will often be detectable and improve the knowledge of 
the lens configuration; but high sampling networks will be needed in the case 
of very fast variations during the caustic crossing of double lens events. 



5.3 Limitation from the knowledge of the source distance distribution 

The Magellanic Clouds are stellar populations located at well defined dis- 
tances. Therefore, the allowed parameter space of an observed lens configu- 
ration is less extended than when the source distance is poorly known. The 
distance distribution of the Galactic bulge population is wider, but the rel- 
ative uncertainty on the bright source's positions is less than 10% [S2]. This 
uncertainty does not have a strong impact on the optical depth estimates. In 
contrast, monitored sources at large Galactic longitude can span a wide range 
of distances (see Sect. 16.21) . The interpretation of the mean optical depths re- 
quires either a good knowledge of the source distance distribution (through 
complementary observations such as spectroscopy) or a model that allows the 
simulation of the survey catalogs and to predict the mean optical depths [57] . 



4 Probing Lensing Anomalies NETwork 

5 Microlensing Follow-Up Network 

6 Microlensing Network for the Detection of Small Terrestrial Exoplanets 

7 Large Synoptic Survey Telescope 
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6 Review of results 

6.1 Galactic model 

We describe below the model and the parameters of one of the simplest Galac- 
tic models, that is commonly used for the optical depth interpretations. To 
facilitate the comparisons, we will interpret all the optical depth results and 
some event duration distributions within this framework, and briefly mention 
the alternatives in the subsequent discussion. The model includes a central 
bulge, a thin disk, a hypothetical thick disk and a halo; since the exact mor- 
phology of the Galactic Spiral Arms (GSA) is yet unknown [49 , no spiral arm 
features have been included in the model. The parameters we use here are 
summarized in Table [2] 

— The observer position and velocity: To calculate the optical depths 
and duration distributions within the Galactic plane, one needs to take 
into account the position of the Sun with respect to the disk (15.5 ± 3 pc 
above the median plane, according to [6 3) ) and the peculiar solar motion 
with respect to the Local Standard of Rest [35] : 

v QR = 10.0, v Q g = 5.3, v Qz = 7.2 (kms- 1 ). (10) 

— The Galactic halo: Many models of the Galactic halo have been tested 
to interpret the observations; from the simple "standard" spherical halo 
to flattened halo. Discussions can be found in [3T] [17J [53], and the corre- 
sponding predicted optical depth maps can be found in |50) . We will only 
consider here the so-called S-model [12] that allows the easiest comparisons 
and combinations between the various published results towards LMC and 
SMC. This spherical halo is considered as isotropic and isothermal with a 
density distribution given in spherical coordinates by : 

R 2 Q + R 2 C 

p H (r)=PhQ r 2 +R2 > ( n ) 

where phQ is the local halo density, Rq is the Galactocentric radius of 
the sun, R c is the halo "core radius" and r is the Galactocentric ra- 
dius. A Maxwellian distribution is assumed for the velocities with a one- 
dimensional dispersion am- Usually authors consider mono-mass distribu- 
tions for the lenses when simulating microlensing. 

— The Galactic disks: The disk densities are modeled by a double expo- 
nential expressed in cylindrical coordinates: 

S e f -(R-R Q ) \ f-\z\\ 
P n(R,z) = —exp^ Jexp(^ — j , (12) 

where Sq is the column density of the disk at the Sun position, H the disk 
height scale and h its length scale. A thick disk has been discussed as a 
possible component of the Galactic structure [S5] [53] [10] [H] [HI] , for the 
interpretation of the optical depth towards LMC/ SMC and far from the 



15 



Table 2 Parameters of the Galactic model. 





Rq (kpc) 


8.5 






Thin 


Thick 




Z Q (M pc- 2 ) 


50. 


35. 




H (kpc) 


0.325 


1.0 




h (kpc) 


3.5 


3.5 


Disks 


M tMn (xlO lo M ) 


4.3 


3.1 




o r (km s^ 1 ) 


34. 


51. 




erg (km S -1 ) 


28. 


38. 




er z (km s _1 ) 


20. 


35. 




a (kpc) 


1.49 






b (kpc) 


0.58 




Bulge 


c (kpc) 


0.40 






Inclination <P 


45° 






M B (xW 10 M e ) 


1.7 






"bulge (km S _1 ) 


110. 






p hQ (M Q pc- a ) 


0.0078 


Halo 


Rc (kpc) 


5.0 




S-model 


exin (km.s^ 1 ) 


155. 






M in 60 kpc (1O 1O M ) 


40 






P Q (MQpC~ i ) 


0.085 


0.098 


Predictions 


Vrot Out sun (km s 1 ) 


211 


222 




TLMC 


4.7 x 10~ 7 






TSMC 


6.58 x 10" 7 





Galactic center. The kinematics of these disks involves a global rotation 
given by [34 



V rot (r) = V ro t.Q x 



1.00762 



Rc. 



0.0394 



0.00712 



(13) 



where Vrot,© = 220 km/s [31], and an anisotropic Gaussian peculiar veloc- 
ity distribution, characterized by the velocity dispersions of Table [2j The 
mass function for the lenses is taken from [57] for both the thin disk and 
the bulge. 

The Galactic bulge: The density distribution for the bulge - a bar-like 
triaxial model - is taken from [17] model G2, given in Cartesian coordinates: 



PB 



6.577ra6c 



-r 2 /2 



© 2 + ©" 



-I 2 



(14) 



where Mb is the bulge mass, and a, b, c the length scale factors. The 
inclination of the bar |101] differs with the authors (from ~ 14° according to 
[55] to 45° according to [55] and 84 ) . The transverse velocity distribution 
of the bulge stars is given by 



Jt{vt) 



-vt exp 



bulge 



^® bulge. 



(15) 
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Fig. 3 V«ei« o/ t/ie Galactic plane and fields monitored by EROS. 



6.2 Results towards the Galactic plane 



Figure [3] shows the fields monitored by EROS towards the Galactic plane. The 
fields monitored by MACHO, OGLE and MOA are all included within the 
bright zones of the Galactic center. MACHO also monitored fields close to 

jScuti {I ~ 18°) and OGLE monitored fields at I 20° and / 30°, but 

only EROS published optical depth measurements far from the Galactic Cen- 
ter > 8°). More than 4000 events have been detected towards the Galactic 
bulge (but only a fraction of them have been used to obtain optical depths un- 
der controlled efficiency), and 27 towards the Galactic Spiral Arms (22 used for 
optical depth determination). The interpretation of the optical depth measure- 
ments is somewhat complicated by the limited knowledge of the lensed source 
distances. Indeed, the optical depth estimated for the full sample of detected 
stars or with DIA analysis is the average over all the stars across the lines of 
sight of the GC. The optical depths towards the spiral arms also correspond to 
average values. The uncertainty on the source distances is reduced towards the 
Galactic center by selecting the clump giant stars, with the triple advantage 
to monitor a well located population (D$ = 8.5 kpc± 10%), less affected by 
blending problems, and benefiting from a better photometric accuracy because 
of the brightness of the sources. The sources monitored towards GSA span a 
wide range of distances (±5 kpc) and their mean distance — estimated to be 
7 ± 1 kpc — is also somewhat uncertain |87j . In the latter case, the measured 
optical depth towards the monitored population corresponds to expression (J7]) 
where dns{Ds)/dDs is the distance distribution of the monitored stars, cor- 
responding to the true distribution distorted by the star detection efficiency. 
Because of this limited knowledge of the source distances, the results are com- 
pared with optical depth estimates at various distances. The following optical 
depth measurements are summarized in Fig. 0] (here < > is corrected for 
the detection efficiencies): 
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events 
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<t E > 


reference 


sons 


deg. 2 


analyzed 
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1° , 5° 


xlO 6 


corrected 


OGLE |112| 


2 


0.81 


all 


9 


±5, -3.5 


3.3 ± 1.2 




MACHO [TTI 


1 


12. 


all 


45 


2.55,3.64 


3 9+ 18 




MACHO [16] 


3 


4. 


all/DIA 


99 


2.68, -3.35 


~\~ ■ 5 2 
OjCO -0.50 




EROS S3 


3 


15. 


bright 


16 


2.5, -4.0 


0.94 ± 0.29 




MOA [1041 


1 


18. 


all/DIA 


28 


4.2, -3.4 


3 36+ 1 11 




MACHO [86] 


7 


4.5 


bright 


62 


1.5, -2.68 


„ +0.47 
^•^'-0.38 


21.6 ± 3 


OGLE \W5\ 


4 


5. 


bright 


32 


1.16, -2.75 


z - oo -0.46 


28.1 ± 4.3 


EROS [62] 


7 


66. 


bright 


120 


GC 


map 


28.3 ± 2.8 


EROS IFH 


7 


20.1 


all 


22 


GSA 


map 


48. ± 9. 



Many of these measurements are correlated because of the large overlap 
between the survey's exposures (N b s x AT b a ). Merging these results would 
require a close collaboration between the authors of these surveys. 

The variation with the latitude deduced from the largest sample (EROS) 
is well fitted by 



t/1(T 6 = (1.62 ±0.23)exp[-a(|6|-3°)], with a = (0.43 ±0.16) deg" 1 . (16) 

This fit agrees with the latest results of MACHO [86] and OGLE-II [TPS] 
that use only bright stars, and with the Galactic models of [50] and [32]. The 
estimates deduced from DIA analysis as well as from searches using all the de- 
tected stars correspond to average optical depth along the line of sight of the 
GC. For a relevant comparison with the optical depth involving bright stars 
only that are all supposed to be close to the Galactic center, we plot in Fig. 
|4] the optical depth Tbuige corrected for the parasitic disk-disk lensing. Nev- 
ertheless, these "all star measurements" seem to be somewhat systematically 
larger than the measurements restricted to bright stars. The corrections for 
disk-disk lensing may be a little too large, or the impact of disk bright stars 
may have been underestimated. 

A local excess of optical depth has been reported by [SB] with a 2cr statisti- 
cal significance, but no conclusive evidence for clustering of events was found. 
The measured tqsa agrees with the predictions of the simplest Galactic model 
(Fig. U]down), providing no evidence for a thick disk. 

Detailed maps of the expected microlensing optical depth are now available 
[69] for the interpretation of the available and future observations towards the 
Galactic bulge. They provide optical depths up to a given magnitude. 

The evaluation of the efficiencies needed by the optical depth determi- 
nations allows us to correct any distribution in order to extract the prior 
distributions. This has been used firstly to check the statistical properties of 
the events as already mentioned in Sect. 14.51 and then to obtain their prior 

distributions. These £e distributions given by the three collaborations have 
been recently analyzed by [37] to constrain the slope of the Galactic Bulge 
mass function (~ m~ a ) in the main sequence (MS) range; the result from the 
MACHO data set {q-ms = 1.7 ± 0.5) was found compatible with the numbers 
deduced from the OGLE and EROS data sets, and also in agreement with 
previous estimates [122] . 
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Fig. 4 (up) The measured optical depth in the Galactic plane as a function of the absolute 
Galactic latitude \b\. The line shows the fit lldj) . (down) The optical depth as a function of 
Galactic longitude. Results at large longitude \l\ > 8° come from the specific EROS GSA 
search at b ~ —2.5°. Results around longitude zero (MACHO: open circles, EROS: filled 
circles, OGLE: crosses) come from the Galactic center searches at \b\ ~ 2.5°. The lines 
show the predicted optical depths as a function of latitude at 6,7 (thick line), 8, 9, 10, 11 
and 12 kpc (from lowest to highest curve) at b = —2.5°, according to the model described 
in Sect. \6.1\ The measured optical depth around I = 0° are compatible with the expected 
optical depth at 8.5 kpc; the best estimated mean distance for the sources in the Galactic 
Spiral Arms is confirmed to be 7 kpc. 



As a conclusion regarding the Galactic plane directions, there is now a 
satisfactory agreement between the optical depths and the efficiency corrected 
duration distributions obtained by the various teams. This agreement, as well 
as the agreement with the Galactic models, give confidence in the reliability 
of the technique, which is necessary to address the more controversial subject 
of the microlensing towards the Magellanic Clouds. 
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Fig. 5 Approximate positions of the monitored LMC fields. The red rectangles (resp. green 
squares) show the OGLE-II (resp. OGLE-III) fields, the blue dashed squares show the MA- 
CHO fields, the red large square is the EROSl-plate field, the small yellow rectangle within 
the LMC-bar indicates the EROS1-CCD field, and the black rectangles show the EROS-2 
fields. 



6.3 Results towards the Magellanic clouds 

Figure [S] shows the fields monitored by the different teams towards LMC. The 
surveys of SMC are much more limited because of the size of this target. 
The optical depths towards these directions have been studied by MACHO, 
EROS and OGLE. Detection of candidates has also been reported by the MOA 
collaboration [80], but without optical depth estimates. 

The main result from the LMC/SMC surveys is that compact objects of 
mass within [10 -7 , 10] x Mq interval are not a major component of the hidden 
Galactic mass (Fig. [BJ). 

Light objects (< 10 _1 M Q ) in the Galactic halo were clearly ruled out as 
early as 1997 by EROS1 [50] which did not find any short time-scale events 
(with tg < 1 day). Heavy objects (1 to 30M©) contribution to the Galactic 
halo was also seriously constrained as soon as 2001 by MACHO [T9"] . 

The contribution of medium mass objects to the halo has been a more con- 
troversial subject, as background sources as well as self-lensing contributions 
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Fig. 6 Published constraints on the fraction f of the standard spherical Galactic halo made 
of massive compact objects as a function of their mass M . The solid line labeled EROS shows 
the EROS upper limit '106], based on the combination of the EROS1 results with the LMC 
EROS2 results. The dotted (incomplete) line takes into account the SMC search direction, 
assuming that the observed candidate belongs to the halo. The OGLE upper limit 11171/ is 
shown in dashed-dotted line. The closed domain is the 95% CL contour for the f value 
claimed by MACHO and the star surrounded by the small rectangle shows the contribution 
of the 2 OGLE events, assuming that they belong to the halo. 



[95] [96] have complicated the data analysis and its interpretation. The MA- 
CHO collaboration has published non-zero halo signal estimates towards the 
densest part of the LMC, but EROS did not find indisputable events within 
its sub-sample of bright stars in its wider LMC search field, and has only 
published upper limits on the halo contribution. To illustrate some of the dif- 
ficulties of the LMC analyses, we have to mention that three of the MACHO 
[SO] and ten of the EROS [ES] LMC/SMC early published candidates have 
been discarded as viable microlensing events because of a significant loss of 
the signal to noise ratio from longer light-curves or because of the detection 
of a second bump (as late as 8 years later [75]). 

6.3.1 Exclusion limits on the Galactic halo 

Four results are mutually compatible: 

— the EROS upper limit [IDS] . 

— the MACHO high-mass upper limit [T9"] . 

— the most recent OGLE-II upper limit [117] . 

— and the first combination of data from different teams [14] (hereafter called 
EROS1-CCD+MACHO) obtained with the EROS1-CCD high sampled 
data (20 minutes) taken in 1991-95 with the 16 CCD camera [SO] GEE] and 
with the MACHO-2year short events search in 1992-94 through a spike 
analysis [TU] . 

The upper limit, considered as a function fs(M), is the maximum contri- 
bution of objects of mass M to the S-model halo optical depth, expressed in 
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fraction /. The maximum contribution fp of a population with a (normalized) 
mass distribution F(m) different than a Dirac function S(m — M), is given by 



We can combine the "zero event" analysis to obtain a stronger constraint 
towards LMC. Including the non-zero event measurements from the SMC or 
from the OGLE-II |117] analysis in such a combination would need a specific 
collaborative work, since the corresponding exclusion limits depend on the 
expected background of non-halo events in each survey (see below). 

Combination of the limits on the low mass side: The simplest way to 
combine the 95% exclusion results is to consider the lower envelope of the 
exclusion curves from all the null analyses. A more sophisticated technique 
would be to calculate the exclusion limits by adding the signals expected by 
each analysis, taking into account only once those objects that arc monitored 
simultaneously by several surveys. Figure [7] (top panel) shows the numbers 
of events expected from the S-model halo lenses for each of the "zero event" 
experiments (EROS1-CCD, EROSl-plates, EROS2, EROS1-CCD+MACHO, 
MACHO high-mass) @. A complete combination, following 14! would need a 
close collaboration between the teams involved. Nevertheless, a simple combi- 
nation can be made: 

— let n ERO si-piates(M), n E ROSi-CCD+MACHo(M) and n ERO s2(M) be the 
expected numbers of events as a function of the lens mass M for the cor- 
responding analysis; 

— since the EROSl-plates data overlap with the MACHO-2year observations, 
they overlap with the combined EROS1-CCD+MACHO analysis. There- 
fore, the total number of expected events from the EROSl-plates plus the 
EROS1-CCD+MACHO analysis is not the sum of the expected numbers 
but it is at least larger than n E ROSi- P lates (M) and uerosi-ccd+macho {M) 
(grey dashed line in Fig. 0-top panel); 

— since the data sets of the EROSl-plates and of the combined EROS1- 
CCD+MACHO analysis do not overlap with the EROS2 analysis the func- 
tion 

npROS2{M) + max[n E ROSl-plates(M),nEROSl-CCD+MACHo(M)} (18) 

represents a minimum for the total number of events expected from all 
these analyses (thick line in Fig. [TJ-top panel) . 

Using this combination allows us to improve the exclusion limit between 10~ 5 
and 10 _4 M Q with respect to the envelope of the EROS/MACHO upper limits 
(Fig. (bottom panel)). 

8 in these experiments, the zero event observed are the very final numbers after rejection 
from complementary observations, as discussed just before this sub-section. 




1 



(17) 
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-8 -6 .4-2 2 

logM = 2log (<t E >/70d) 



Fig. 7 The top panel shows the number of expected events towards LMC as a function of 
the lens mass M for the S-model. The EROSl-total line shows the sum of the expectations 
from the CCD-camera and from the plates (dashed lines). The MACHO-high mass line 
corresponds to the "zero event" high mass search. The EROS1-CCD+MACHO line (blue) 
is the combined result from the EROS1 CCD-camera and the MACHO 2 year analysis. The 
grey dashed line is the envelope of the EROSl-total and the EROS1-CCD+MACHO lines. 
The final combination (thick line) is the sum of EROS2, of the grey dashed line and of half 
of the MACHO-high mass expectations (see text). In the lower panel, the thick line shows 
the combined 95% CL upper limit on f = tlmc I ^--7 x 10~ 7 based on no observed events. 
The thin line is the EROS1+EROS2 LMC limit. 



Combination of the limits on the high mass side: Another combination 
can be made on the high mass side with the MACHO high-mass "zero event" 
analysis [19] and the EROS results |106] ; the MACHO search did not overlap 
with the EROS searches for more than 3 years (1994-96), in the middle of the 
MACHO 5.7 year analysed observations; therefore one can conservatively add 
half of the total published expectations from the MACHO high-mass search 
to the expectations from the other "zero-signal" searches. This is conservative 
as the detection efficiency — specifically for long duration microlensing — is 
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better for events that have their maximum in the middle of the observation 
period. 

The OGLE upper limit: No simple combination can be made with the 
OGLE-II results |117j . as their field and observation epoch are entirely covered 
by the EROS2 program. As for MACHO, the OGLE-II catalog of monitored 
stars contains fainter stars than the EROS2 catalog, and a combination of 
these results could potentially improve the limits but would need a close col- 
laboration between the teams. Such a combination would also be much more 
complicated, as the OGLE analysis is not a "zero-signal" analysis; like the 
one SMC-event found by EROS and MACHO [2] [13], the 2 events found by 
OGLE-II cannot be unambiguously affected to self-lensing. Because of the ex- 
istence of these events, according to standard statistics, the OGLE-II limit 
should not be estimated by considering that 3 events correspond to the 95% 
CL excluded mean value; indeed, the 95% CL upper limit on the mean signal 
S contaminated by a background known by its expected mean value B, were 
S + B follows a Poisson law, is not 3 when the observed number is equal to B. 
In the present case, if one assumes that B — 2 (which is probably a maximum 
according to [TT7] [§T] Q7] [H] and UM), then following [52], the 95% upper 
limit when 2 events are observed is S < 4.72. This limit would be lower only 
if B > 2. If B = 2, the OGLE upper limit on the halo component should be 
corrected by a factor 4.72/3 = 1.57, but this needs to be finalized using the 
exact value of B. The statistical procedure described in [20] could also be used 
for the OGLE-II non-zero event analysis. 

6.4 The Galactic halo 

6.4-1 The apparent discrepancies towards the LMC 

We can quantify the incompatibility between the MACHO results from [T2] 
and [17] and the EROS results pH]. The 13 observed events by MACHO in [17] 
correspond to a fraction /macho = 0.21 of the S-model halo; The MACHO 
collaboration has estimated that at least 9.6 of these events should be due 
to hidden baryonic matter (Milky- Way or LMC halo). This corresponds to 
9.6/ f macho — 45.7 events expected with a full S-model halo entirely made of 
objects of mass O.6M . For the same mass, about 39 x j "macho = 8.2 events 
are expected in the EROS1+EROS2 LMC data while zero were observed. 
For a mass 0.9M Q , which is la larger than the central MACHO value, 30 x 
/macho = 6.3 events were expected in the EROS data. Assuming an optical 
depth dominated by the Galactic halo and no correlation between the data 
sets, the EROS and MACHO measurements are statistically incompatible. 
The correlations due to the (actually small) overlap between the data sets 
makes this incompatibility worse. 

There are considerable differences between the EROS, MACHO and OGLE 
data sets that may — at least partially — solve the apparent conflict. Generally 
speaking, MACHO uses faint stars in dense fields (1.2 x 10 7 stars over 14 deg. 2 ) 
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while ER0S2 uses bright stars in sparse fields (0.7 x 10 7 stars over 90 deg. ). 
As a consequence of the use of faint stars, only two of the 17 MACHO (set 
B) candidates were sufficiently bright to be compared to the EROS bright 
sample; as the corresponding events occurred before the EROS data taking, 
there is no evidence for a malfunctionning experimental setup during one of 
the surveys. The use of dense fields by the MACHO group also suggests that 
the higher MACHO optical depth may be due, in part, to self-lensing in the 
inner parts of the LMC [61]. The contamination of irregular variable objects 
faking microlensing in low photometric resolution events could also be stronger 
in the sample of faint stars used by MACHO. As already mentioned, another 
difficulty when using faint source stars originates in the blending effects that 
complicate the optical depth calculations. 

6.4.2 The SMC case 

In the framework of the S-model halo, the EROS upper limit on the halo 
fraction from the combined LMC/SMC studies in [4 and [106 assumes that 
TSMC-haio = 1-4 x TLMC-halo'i this would be larger with a flattened halo 
or a thick disk. The dotted line in Fig. [6] shows the small impact of the SMC 
measurements according to this model. Even though the impact on the optical 
depth estimates seems limited, the compatibility of the SMC microlensing 
events population with the LMC one is the subject of an interesting debate. 
The long durations of the four SMC events published in the 5-year search of 
EROS [4, is found to be incompatible with the £e distribution of the MACHO 
LMC events [17], and make these events more compatible with unidentified 
variable-type stars or with SMC-self-lensing. This issue is discussed by [9"T] . 

Because of their long duration, a fair fraction of the SMC events could 
benefit from complementary studies, thus providing indications on the lens 
population. The first event towards SMC (MACHO-SMC-97-1/EROS-SMC1) 
was discovered in 1997 and observed by both MACHO [13] and EROS [55] , 
Thanks to its long duration (£e = 123 days) and its brightness, a parallax 
analysis was performed on the light-curve, and from the absence of measurable 
distorsion, it was concluded that the lens probably belonged to the SMC itself. 
This interpretation was confirmed from the Spectroscopy of the lensed star [57] . 
The second event discovered towards SMC (MACHO-SMC-98-1) was the first 
event from the Magellanic Clouds to benefit from an intensive international 
collaboration between MACHO [T5], EROS BQ, PLANET [S], OGLE [TH] and 
MPfJ! [94 collaborations on 18 June 1998, as it was early detected by MACHO 
as a binary-lens microlensing event. The common conclusion from the different 
survey teams was that the proper motion of the lens was consistent with the 
lens beeing within the SMC rather than in the Galactic plane or halo. But at 
least one SMC-event, OGLE-2005-SMC-001, measured both from Earth and 
from space [J5] , had a parallax estimate such that the probability for the lens 
to belong to the halo was 20 times larger than to belong to the SMC. The two 
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later events do not enter in any optical depth calculation, and they are only 
mentioned here to illustrate the problem of the lens location. Unfortunately, 
the very small statistics limits the interpretations, but we can conclude that 
the SMC lens population is probably not dominated by the Galactic halo, and 
is certainly different from the LMC lens population. 

6.4-3 A synthesis attempt: halo versus local structures 

Taking all the observational results at their face value, we may conclude that 
the initial hypothesis of an optical depth dominated by the Galactic halo — 
almost uniform through all the monitored LMC fields and 1.4 times larger 
towards SMC — is wrong, because it cannot explain the EROS-MACHO dif- 
ferences nor the LMC-SMC differences. 

As the results found by the teams are mutually compatible towards the 
Galactic plane directions (Sect. 16. 2p . the falseness of this hypothesis should 
now be favored rather than the eventuality of wrong measurements. 

Abandoning the original microlensing survey paradigm of a dominant halo 
allows one to explain the variety of the observations by the large differences in 
the monitored fields. Fig. [5] supports this view. If — as suggested by their small 
contribution to the halo — one assumes that lenses do not trace the Galactic 
halo hidden matter, then the apparent discrepancy between the surveys can 
be understood by considering local structures inducing self-lensing; such struc- 
tures (populations of foreground lenses as well as populations of background 
sources [18) 1 may be responsible for the variability of the mean optical depth 
with the monitored fields. Nevertheless, a convincing proof of the existence of 
such structures is still to be provided. One of the aims of the current and next 
generations of microlensing surveys towards LMC (like SupcrMACHO) is to 
improve the knowledge of the lens populations, and specifically their location. 

For a pertinent synthesis, it seems that the total optical depth — understood 
as an average over the monitored fields — is the perennial observable that will 
allow the building of contour maps of the massive compact object spatial dis- 
tribution. The following numbers have to be imperatively associated with the 
corresponding fields: 

— tlmc-macho — (1.0±0.3) x 10~ 7 , after revision by [3U]- There have been 
some controversy on the MACHO optical depth estimates. The first esti- 
mate from [T2] was clearly much too high (r = 2.9^q 9 x 10~ 7 ). This result 
was seriously downsized to r = 1.2^ ^ x 10 -7 after the analysis of 5.7 years 
of data [T7]. MACHO collaborators and outsiders have revised/disputed 
the later result; |30| took into account known contamination of the signal 
by non-microlensing variabilities and published a revised optical depth of 
r = (1.0 ± 0.3) x 10~ 7 . [28 and [29] have re-analized the publicly avail- 
able MACHO light-curves with a neural network system and estimated a 
number of contaminants larger than the one quoted in the original MA- 
CHO paper. They conclude that Tlmc-magho — (0-4 ± 0.1) x 10~ 7 . A 
controversy followed ([BIT] and [ST]). 
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- tlmc-eros < 0.36 x lO" 7 (95% CL). 

— tlmc—OGLE = (0.43 ± 0.33) x 10~ 7 is the total measured optical depth 
(taking into account 2 observed events). 

The distributions of event durations (except for EROS which has no event) 
have also been used by several authors to refine the interpretations. For ex- 
ample, a detailed study of the consequences of the MACHO results was per- 
formed in [68 , which takes into account the optical depths, rates and durations 
expected from several Galactic structures (thick disk, spheroid) through the 
mass moment method [32] [35] , and concluded that at least 3 of the observed 
microlensing events were due to LMC lenses, but that the others were prob- 
ably due to other Galactic structures. More recently [36] concluded that up 
to about half of the MACHO events could be attributed to the LMC dark 
matter halo instead of the Milky- Way dark halo. Calchi Novati et al. [55J also 
conclude from their rate and durations that the 2 events reported by OGLE- 
II are perfectly compatible with the expectation from known lenses {i.e. not 
dark), dominated by the LMC self-lensing. 
Conclusion 

A substantial contribution of compact objects to a standard halo is now 
clearly excluded. Compact objects belonging to a flattened halo or a thick 
disk are also seriously limited by observations both towards SMC and GSA. 
Wether or not there is a small dark matter component under the form of 
compact objects is still an open question. 



7 The future of the optical depth studies 

Several authors have investigated the possibility to systematically use comple- 
mentary observations — instead of occasionnally as currently done — to look 
for non-standard features. Systematic searches for parallax effects towards 
LMC have been examined by [88] based on the use of the ESO-NTT tele- 
scope, and by |103j using measurements by the ESO-VLT or by the Hubble 
Space Telescope on strongly magnified sources. Such follow up measurements 
would increase the information about the events and help to localize the lens 
populations. Accurate mapping of the optical depth is another way to improve 
the knowledge on these populations and to measure the self-lensing contribu- 
tions. Microlensing searches that address this question are currently running 
or about to observe: SuperMACHO [35], OGLE-III and IV [TUT] and the in- 
frared VVV project |114j using VISTA facility, should renew this scientific 
domain within the next few years. The LSST program will certainly dominate 
the following decade of microlensing studies as a microlensing factory, possibly 
with contributions from spatial missions like GAIA. 
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